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Abstract. In this paper we propose an algorithm that correctly discards a 
set of numbers (from a previously defined sieve) with an interval of integers. 
Leopoldo's Theorem states that the remaining integer numbers will generate 
and count the complete list of primes of absolute value greater than 3 in the 
interval of interest. This algorithm avoids the problem of generating large lists 
of numbers, and can be used to compute (even in parallel) the prime counting 
function tt (h). 
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1. Introduction 

An old problem in mathematics is the way to compute the amount of prime 
numbers less or equal to a given value h [21 p. 347]. This function is known as tt (h) 
[4]. The preeminent method for such task since the 3rd century BC was the sieve 
of Eratosthenes. From then on there were no great advances on the subject until 
Gauss in 1863 [21 P- 352]. His work allowed more advances 0|6] and then boosted 
by the growth of calculation power in the 20th century. Recent implementations 
[6] require several sophisticated computational strategies. 

One of the major difficulties in solving this issue -or may be the only one- is 
the question of how prime numbers are distributed among integers. Matheaticians 
and other scientists have carried great efforts to find some general formula and 
simple prime generating formula. Such formula is still ellusive. Some formulas have 
been proposed, but they are valid only in certain cases, such as Mersenne numbers, 
and others. The key to find a simple formula that would allow us determine all 
prime numbers, should begin demonstrating that all numbers obey a general logic 
order, from there on it could be possible to find a sequence that would permit us 
to establish such formula. In other words it is about establishing a new sieve of 
universal validity. This basic essential aspect was the one which guided L. Euler in 
1770 to express "Mathematicians have tried in vain to discover some order in the 
sequence of prime numbers but we have every reason to believe that there are some 
mysteries which the human mind will never penetrate". 

i 
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On the other side, the efforts of mathematicians like Mersenne, Fermat, and 
many others from Ancient Greece up to our days, to introduce certain logic in the 
calculus of some prime numbers with simple arithmetic operations became useless 
due to the observation made by Gauss [2) himself, who stated that the amount of 
prime numbers less than h is near the value of the integral: 



This observation allowed Gauss to introduce the function tt (h) toghether with a 
probabilistic argument in the studies about the distribution of primes. In fact the 
function tt (h) is related with the Fehlerintegral or Error function (Erf) , introduced 
by Gauss to analyze the uncertainties that appear in experimental measurements. 
It is to be noted, that in experimental sciences like Physics or Biology, the scientist 
needs to justify that the results of a certain experiment can be treated inside the 
gaussian frame. This a priori it does not seem to happen in the mathematics area, in 
particular in Number Theory, were Axioms play a fundamental role. Nevertheless, 
in the last years a lot of works have merged regarding the issue of the location of 
prime numbers, whether and if primes are located in a random way or if they have 
a totally determined structure. The papers that were published in this area are 
various and interesting (see references [9], [TU1 [TTJ [131 [TBI ■ These results are not 
completely satisfactory, therefore "the human mind" -according to L. Euler- will 
continue the attempt to penetrate the mysteries of the apparent unbreakable tower 
of primes. In particular, about this matter, Terence Tao [8] proposes the existence 
of certain dichotomy between structure and randomness among primes, and tries to 
separate the structured component (or correlated) from the "pseudorandom" com- 
ponent (or decorrelated) . Anyway, the existence of a distribution of prime numbers 
affected by the dichotomy of structured-randomness implies that the structured 
behavior or low complexity behavior can be determined by any general formula, 
meanwhile in the high complexity behavior we could recognize at least two aspects: 
i) A totally randomness aspect or decorrelated, and the other, ii) A chaotic one, 
which obeys some kind of logic formulation. In our opinion any well known prime 
number allows us to describe a sieve, which can be enlarged as far as new primes 
are verified as such numbers [Q. However we consider that the unknown prime 
numbers, are not undefined nor undetermined. That means that all of them can 
be defined and determined by a general and unique algorithm, logically determined 
and located sequentially among integers [5]. 

Returning to our primary purpose of implementing an algorithm to calculate the 
value of the 7r(h) function for a given h, in this paper we will study this problem 
defining a new sieve whose properties permit an elementary study of ir(h), toghether 
with the possibility of finding its value on a given interval. The basis of this sieve 
stands on an infinite matrix presented in [I], which has the property of generating 
all indexes that lead to non prime numbers according to the form 6n + 1 or 6n — 1 , 
a similar approach can be found in [7], but with a different implementation. 

The paper is organized as follow. In Section 2. we reviewed the properties of 
the form 6n + 1, being n = 0, 1, 2, . . . , which is the base of our research. Section 
3 is devoted to our proposal of a sieve to define and compute all prime numbers 
larger than 3. In Sections 4 and 5 we present the Leopoldo's Theorem and the 
possibility to compute Tt(h) function according with the A algorithm that we are 
introducing. Section 6 and 7 are devoted to our actual results, comments on future 
research work and conclusions. 
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2. About the form 6n + 1 

In [T] we reviewed some well known properties of numbers: 

Theorem 2.1. Every prime number of absolute value greater than 3 can be written 
in the form 6n + 1 or 6n — 1 . 

Proof. Let's see the equivalences modulo 6. Suppose g prime. 

1) If g = =*>g = 6n^6 | q, ABS. 

6 

2) If g = 1 =>g = 6n+ 1, which is not impossible since 7 = 6 + 1 is a prime. 

3) If q = 2 =*>g = 6n + 2 = 2 (3n + 1)=^2 | g, ABS. 

4) If g = 3 =^g = 6?i + 3 = 3 (2n + 1)^3 | g, ABS. 

5) If g = 4 =^g = 6n + 4 = 2 (3n + 2)=^2 | g, ABS. 

6 

6) If g = 5 =>g = 6?i + 5 = 6n + 6 — 1 = 6 (n+ 1) — 1, which is not impossible 

6 

since with n = 1 this gives 11, a prime. 

Then we introduce some definitions: □ 

Definition 2.2. The a class of integer numbers [5] is the set 
(2.1) a = {xeZ/x = 6n + l,neZ} 



Definition 2.3. The (5 class of integer numbers [5] is the set 

(2.2) /3 = {xeZ/x = 6n-l,neZ} 

The different values of relations in (|2.ip and l|2.2p are shown in Table [TJ 



n 


/3 n = 6n - 1 


a n = 6n + 1 








-5 


-31 * 


-29 * 


-4 


-25 


-23* 


-3 


-19 * 


-17 * 


-2 


-13 * 


-11 * 


-1 


-7 * 


-5 * 





-1 


1 


1 


5 * 


7 * 


2 


11 * 


13 * 


3 


17 * 


19 * 


4 


23 * 


25 


5 


29 * 


31 * 









Table 1. Values of 6n — 1 and 6n + 1 for several n. With (*) we 
mark prime numbers. 



Strictly speaking, a "complete" list of all prime numbers of absolute value greater 
than 3 ({•■■, —7, —5, 5,7,.. .}) is the list of primes from both classes, a and /3. 
We now state a property given in [5] and [I]: 
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Theorem 2.4. Every prime number of absolute value greater than 3 (except for 
the sign) is generated by 6n + 1, with n integer. 

Proof. We must prove the equivalence (except for the sign) between both families 
given in Theorem 12.11 

Let be f a (n) = 6n + 1 and fp (n) = 6n — 1, we must now prove that f a (—n) = 
—fp (n). Indeed: 

f a (-n) = 6 (-») + 1 = -6n + 1 = - (6n - 1) = -fp (n) 

□ 

Definition 2.5. We define the set of integer numbers G a : 

(2.3) G a = {g e Z/6g + 1 is a prime} 

This means, G a is the set of all numbers that (except for the sign) generate all 
primes of absolute value greater than 3 by the relationship (|2.ip . 
In other words, this result is perfectly logic because it is based on: 

i) Peano 's Axioma applied to both sequences of numbers a and (3 starting 
in 1 and 5, respectively, taking into account that the next number in 
those sequences are defined as 1 + 6, and 5 + 6, respectively, and so on 
and so forth; and 

ii) Remark (|3.5p which follows from Theorem 3.4., represented by the Ta- 
ble of Products [4l This table states that products between numbers 
belonging to the a sequence and products between numbers belonging 
to the (3 sequence, results in numbers of Class a; while the crossed prod- 
ucts between numbers of the sequence a with numbers of the sequence 
(3, results in numbers of Class (3. Besides, note that crossed products 
between numbers belonging to the six classes we presented in Theorem 
2.1. never result in numbers belonging to classes a or j3. Table [2] shows 
the results of these products from n — to n = 10; remember, prime 
numbers greater that 3 are only allocated in classes a or (3. 



n 


6n + 1(a) 


6n + 2 


6n + 3 


6n + 4 


6n + 5 {fi) 


6n + 6 





1 


2 


3 


4 = 2x2 


5 


6 = 2x3 


1 


7 


8 = 2x4 


9 = 3x3 


10 = 2 x 5 


11 


12 = 2 x 6 


2 


13 


14 = 2 x 7 


15 = 3 x 5 


16 = 2 x 8 


17 


18 = 2 x 9 


3 


19 


20 = 2 x 10 


21 = 3 x 7 


22 = 2 x 11 


23 


24 = 2 x 12 


4 


25 = 5 x 5 


26 = 2 x 13 


27 = 3 x 9 


28 = 2 x 14 


29 


30 = 2 x 15 


5 


31 


32 = 2 x 16 


33 = 3 x 11 


34 = 2 x 17 


35 = 5 x 7 


36 = 2 x 18 


6 


37 


38 = 2 x 19 


39 = 3 x 13 


40 = 2 x 20 


41 


42 = 2 x 21 


7 


43 


44 = 2 x 22 


45 = 3 x 15 


46 = 2 x 23 


47 


48 = 2 x 24 


8 


49 = 7 x 7 


50 = 2 x 25 


51 = 3 x 17 


52 = 2 x 26 


53 


54 = 2 x 27 


9 


55 = 5 x 11 


56 = 2 x 28 


57 = 3 x 19 


58 = 2 x 29 


59 


60 = 2 x 30 


10 


61 


62 = 2 x 31 


63 = 3 x 21 


64 = 2 x 32 


65 = 5 x 13 


66 = 2 x 33 

















Table 2. Numbers belongning to the six classes we divided nat- 
ural numbers from n = to n — 10. Prime numbers are shown in 
bold letters. 



Definition 2.6. 
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3. A PROPOSED NEW SIEVE 

Definition 3.1. Let A be an infinite matrix whose element a(i,j) Qis 
(3.1) a(i,j)=i + j(6i + l) 

where i, j £ Z. 

In Table [3] we show the elements in the central part of the A matrix. Note that 
axis numbers also match this representation. 



-96 


-71 


-46 


-21 


4 


29 


54 


79 


104 


-73 


-54 


-35 


-16 


3 


22 


41 


60 




-50 


-37 


-24 


-11 


2 


15 


28 






-27 


-20 


-13 


-6 


1 


8 








-4 


-3 


-2 


-1 





1 


2 


3 


4 


19 


14 


9 


4 


-1 










42 


31 


20 




-2 










65 

88 


48 






-3 
-4 











Table 3. Elements of the central portion of A 



3.1. Properties. 

Theorem 3.2. A is symmetrical. 

Proof. A simple expansion shows that 

a (i, j) = i + j (6i + 1) = i + 6r/ + j = j + i (6j + 1) = a (j, i) 

□ 

Definition 3.3. Let A be the set of unrepeated elements of A excluding the axes 
(the elements of the form a (i, 0) y a (0, j)). 

3.1.1. About the signs of the elements of A. Four quadrants can be distinguished 
as shown in Figure [XT] 

II t I 

iv | III 
Figure 3.1. Quadrants of A 

What happens to the signs of the elements of A from each quadrant? 
Because of Theorem 13. 2^ we should only focus on the sign of elements of A 
originally from quadrants I, II, and IV. 

• In quadrant I (i > 1, j > 1) all elements are positive 

• In quadrant II (i < —1, j > 1) 

It's easy to see that j (6i + 1) < 0, then a < Wi, j. 



Coordinates are in the Cartesian sense. 
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• In quadrant IV (i < — 1, j < — 1) 

i + j (6i + 1) = i + j + 6ij. (i + j) < — 1 y ij > \i + j\ > 1. Then, the sign 
is positive. 

Theorem 3.4. The elements of A DO NOT generate prime numbers. 

Proof, a = i + j (6i + 1) with i ^ and j ^ 0. If we put this into (|2.ip and 
suppose p prime, then 

p = 6a(i,j)+l =6(i+j (6i + 1))+1 = 6i+6j (6i + 1)+1 = (6i + l)+6j(6i+ 1) = (6i + 1) (6j + 1) 
but since i and j are different from zero, then p would be a composite, ABS. □ 

Remark 3.5. According to the signs of i and j, 6a + 1 sweeps (except for the 
sign) all possibilities: 

(1) If i > y j > 0, the generated number is of the form a ■ a. 

(2) If i < y j < 0, the generated number is of the form [3 ■ (3. 

(3) If i y j have opposite signs, the generated number is of the form a ■ (3. 
See Table |4] for properties of products of a's and /3's. 



X 


a 


a 


a 


a 


a 




a 


a 



Table 4. Table of products. 



Then, summarizing previous results, we can state that for positive prime num- 
bers: 

(1) The only even prime number 2 belongs to the sequence defined by the form 

2 + 6n for n = 0. All the other even numbers in the same sequence are 
multiples of 2, then they are composites. 

(2) All the others even numbers that belong to sequences of the form 4 + 6n 
and 6 + 6n, being n = 0, 1, 2, are multiples of 2 and obviously they are 
composites. 

(3) The first odd prime number 3 belongs to the sequence defined by the form 

3 + 6n for n — 0. All the other odd numbers in the same sequence are 
multiples of 3, then they are composites. 

(4) Then, except 2 and 3, all the other prime numbers belong to classes a 
(l + 6n) or (3 (5 + 6n). 

(5) Therefore, taking into account Table of Products a x a = a, (3 x (3, and 
a X (3 = j3, all composite numbers generated by these products belong to 
classes a and (3. 

(6) Then all a prime numbers are those expressed by (l + 6n) for n = 0, 1, 2, 
except: 

(a) Number 1 expressed by 1 + 6n, being n = 0. 

(b) All even and odd powers of all class a numbers, as (1 + 6n) 2+k , for 
n = 1,2,3, and k = 0,1,2, ... 

(c) All products between (3 numbers and all their even and odd powers, 
as: [(5 + 6n)(5 + 6(n + m))] 1 +' £ , for n = 0,1,2,..., m = 0,1,2,..., and 
ft = 0,1, 2,... 

(7) All (3 prime numbers, except all the products of each (3 number with all 
a numbers, except 1, like (5 + 6ri)(l + 6m), for n = 0, 1, 2, and m = 
1,2,3,4,... 
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4. Leopoldo's Theorem and the tt function 

4.1. Leopoldo's Theorem. We defined the A set as a list of all the non repeated 
off-axis elements of A. A simple expansion sowed that the elements of A do not 
generate prime numbers by l|2.ip . Finally, we stated and proved Leopoldo's theorem 
ED: 

Theorem 4.1. (Leopoldo's TheorerM) G a =Z- A. 

Proof. Every prime number of absolute value greater than 3 is either a or f3, and 
products between those classes of equivalence are closed on themselves. Because 
of Theorem 13.41 and Remark 13. 5[ we know that A generates all possible a and j3 
composite numbers, and so the elements n of Z ^ A generate prime numbers (except 
for the sing) by 6n + 1. □ 

This means that all integers not generated by l|3.ip will generate all primes of 
absolute value greater than 3 by l|2.ip (and thus, this works as a sieve). 

In order to see the whole picture we are going to suppose that instead of integers 
we are using real numbers for i and j. This will lead us to a graphical sight of the 
matrix where, keeping a(i,j) as an integer, we should see hyperbolas. The plots 
will represent any integer c generated by the matrix, and they can be considered 
as level sets. These sets have certain coincidences when i and j are integers at the 
same time, obviously this is a point in which the elements in the matrix coincide 
(in space) with an element in the set of curves. Some of the curves have lots of 
repetitions specially when c is large, that is why the calculation of all the elements 
of A is futile. What we are looking here, in order to find primes, are the integers c 
which have integer generators i and j. This will lead us to composite numbers and 
indirectly to primes. Some level sets of A are shown in Figure 1570 

Now, suppose you wish to calculate all prime numbers of absolute value greater 
than 3 up to a certain value h = 6c + 1 (c > 0), this means computing tt (h) using 
Leopoldo's Theorem. At first sight, one would have to: 

(1) Generate all elements of A up to c which means that \a = (h — 1) /6 

(2) Discard the axis elements 

(3) Sort the rest 

(4) Discard repetitions 

(5) Remove them from the interval [— c, c] 

(6) Count the remaining numbers 

(7) Apply l|2.ip to show the primes in the interval 

With large numbers, this computation would quickly become time and memory 
prohibitive. However, a closer look at the distribution of elements in the sieve, 
gives an appropriate answer to this problem. 

5. Leopoldo's Theorem and the computation of the tt function 

This method is based on discarding non prime generators. In fact we only need 
to find two integers i and j to verify that the number c that we are checking belongs 
to A. Once we find these numbers there is no need for further search. This will 
simplify the problem by avoiding repetitions. In this section we are going to present 
an approach for an algorithmic method to compute it based in the ideas presented 
above. 



2 We named theorem 4.1 as Leopoldo's theorem because it was stated by Leopoldo Garavaglia in 
the Summer 2006, Spain, but the basis was not published until 2007 [5] and it was proved in [T] 
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Figure 5.1. Several level sets of / (x, y) — x + y (6x + 1). 




In Figure l5~2l we show the representation of the level sets / (x, y) — ±c (c > 0). 
We must only consider non-repeated elements of A originally from within the "star" 
delimited by 

(5.1) f(x,y) = ±c = ±(h-l)/6 
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5.1. An algorithmic approach. An exploration of the level sets allows algorith- 
mic approach to tt (h). 

Algorithm 5.1. We define the A algorithm of arguments c\ and c 2 (c% > c\ > 8), 
as the procedure that 

(1) Declares a natural variable L = 

(2) For c taking every integer value from c\ to c 2 

(a) For integers x from x = — |_(e+l) /5J to x = — + 6c + l) /6j 
verifies if 

c — x 
6x + 1 

takes an integer valu^l. 

(i) if it finds one, adds 1 to L and goes to step (2)(c) 

(ii) if it doesn't find any, goes to step (2)(b) 

(b) For integers x from x = 1 to x = + 6c — l) /6j verifies if 

c — x 
6x + l 

takes an integer value. 

(i) if it finds one, adds 1 to L and goes to step (2)(c) 

(ii) if it doesn't, prints c and 6c + 1, and goes to step (2)(c) 

(c) For integers from x — — [(c+ 1) /7\ to x = — 1 verifies if 

— c — x 
Qx + l 

takes integer values. 

(i) if it finds one, adds 1 to L and goes to the next value of c. 

(ii) if it doesn't find any, prints — c and —6c + 1, and goes to the 
next value of c. 

(3) Once the process has been completed up to C2, reports the accumulated 
value of L: 

A (Ci, C 2 ) = L fi na i 

This algorithm tells the amount of numbers that do not generate prime numbers 
by 6c + 1 in the interval [cijCa], and also gives the list of the missing values as 
well as the primes generated by them. So, the amount of prime numbers between 
hi = 6ci — 1 and /12 = 6c 2 + 1 (ci > 8) is: 

(5.2) Att = 2 (c 2 - cO - A( Cl ,c 2 ) + l 
In the particular case of c\ = 8 (hi = 47) and h = 6c 2 + 1: 

(5.3) 7r(/ l ) = 2c 2 -A(8,c 2 ) 

6. Results and Future Works 

In Table [21 it may be seen that the the proposed algorithm for the n function 
agrees with known results for several testing values. Equations l|5.2p and (|5.3p 
enable a parallelization of the computation of the tt function with the A algorithm. 

Calculation time grows with the parameter c. The last value in the table took 
nearly an hour to be computed with an AMD Athlon 64 X2 Dual Core 4200+ . 
All values were calculated using (|5.3[) . The memory used remains stable, since it 
doesn't take more than it needs to store the variables and operations in steps (2) (a), 
(2)(b) and (2)(c). 



3 This means, (c — x) = (modulo 6ic + 1). 
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h 


ir(h) 
(Mathematica 5) 


tt(/i) 
(A Algorithm) 


tt {h - 3) 
L4J 


10 2 + 3 


27 


27 


25 


10 a + 3 


168 


168 


168 


10 4 + 3 


1229 


1229 


1229 


10 b + 3 


9593 


9593 


9592 


lO^ + S 


78499 


78499 


78498 


10 Y + 3 


664579 


664579 


664579 



Table 5. A algorithm based results versus several known values of n (h). 



Given the fact that this algorithm evaluates all values between c\ and C2, as 
well as the independence of the obtained value An with other intervals, a parallel 
expression of the Algorithm is possible using n independent processors where the 
range of values c for each acting processor is 

c 1P = Cl + (^) B , 
ftp = CV + (^-) P„ 

and so on, where #P is the total number of processors and P n is the number of a 
given processor in the cluster. 

Since the computing load is equally divided between computing nodes, the effi- 
ciency of the process asymptotically approaches 1 (each computing node uses nearly 
100% of itself to calculate). The speedup tends to #P as new computing nodes 
are added whenever (c2 — ci) 3> #P. In a future paper we will deal with the case 

(C 2 - Cl) < #P 

7. Conclusions 

We have seen that using a simple classification of numbers we can group primes 
in families. We have also seen that primes and its multiples are grouped in the 
same classes according to the table of products [H This classification leads us to 
the separation of composite numbers among primes, and the characteristic index 
number that do not generate primes. The matrix A which is the matrix of indexes 
(without the axis) and we use it in order to find primes. Also by analyzing the 
matrix as continuous hyperbolas we have found a way to count the amount of 
elements of A from a set of integers. In this way by counting the amount of "non 
prime generating indexes" we can obtain indirectly the amount of primes in a given 
section. The advantage of this method is that we can calculate the amount of 
primes by section. That means that this is not only an accumulative method, we 
can use it in any section of natural numbers and also use various computers to 
process different sections at the same time. Finally we can also use known results 
of 7r(/i) and start counting primes from the last well known result. As mentioned 
above, in future works we would like to analyze the use of several computers in 
order to find ir(h) and test the proposed algorithm. 
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